Quantum renormalization group of XYZ model in a transverse magnetic field 



o 
o 

(N 



A. Langari 

Institute for Advanced Studies in Basic Sciences, Zanjan 45195-159, Iran 
Max- Planck- Institut fur Physik komplexer Systeme, Nothmtzer Str.38, 01187 Dresden, Germany 

We have studied the zero temperature phase diagram of XYZ model in the presence of transverse 
magnetic field. We show that small anisotropy (0 < A < 1) is not relevant to change the uni- 
versality class. The phase diagram consists of two antiferromagnetic ordering and a paramagnetic 
phases. We have obtained the critical exponents, fixed points and running of coupling constants by 
implementing the standard quantum renormalization group. The continuous phase transition from 
antiferromagnetic (spin-flop) phase to a paramagnetic one is in the universality class of Ising model 
in transverse field. Numerical exact diagonalization has been done to justify our results. We have 
also addressed on the application of our findings to the recent experiments on CS2C0CI4,. 
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Systems near criticality are usually characterised by 
fluctuations over many length scales. At the critical point 
itself, fluctuations exist over all scales. At moderate tem- 
peratures quantum fluctuations are usually suppressed 
compared with the thermal ones. However if tempera- 
ture is near zero, quantum fluctuations especially in the 
low-lying states dominate thermal ones and strongly in- 
fluence the critical behavior of system. Zero-temperature 
(quantum) phase transition may occur in the area of spin 
systems by applying noncommuting magnetic field which 
introduces quantum fluctuations. Such a situation has 
been studied in the three dimensional Ising ferromag- 
net L1H0F4 in a transverse magnetic field — . However 
due to its high dimensionality, the system behaves in a 
mean-ficld-like manner. In this paper we are going to 
consider the one-dimensional XYZ model in the presence 
of a transverse field where quantum fluctuations of sym- 
metry breaking field play an essential role. Generally 
Renormalization Group (RG) is the proper method to 
give us the universal behavior at long wave lengths where 
other methods fail to work accurately. 

The spin-(s =)~ Hamiltonian of this model on a peri- 
odic chain of N-sites is 



H 



N 

E 

i=l 



[J x *?<T? +1 + Jya?aV +1 + Aa!a* +1 -ha?], (1) 



where J x > and J y > are exchange couplings in 
the XY easy plane, < A < I is the anisotropy in Z 
direction which is in J y units and h is proportional to 
the transverse field. a a ; a — x, y, z are Pauli matrices. 

When h = 0, the XXZ model (J x = J y ) is known to 
be solvable and critical (gapless) while — 1 < A < l 2 . 
The Ising regime is A > 1 and A < —1 is the ferro- 
magnetic case. Magnetic field in the anisotropy direction 
commutes with the Hamiltonian (h — 0) and extends 
the gapless region (quasi long range order) to a border 
where a transition to paramagnetic phase takes place. 
The model is still integrable and can be explained by a 
conformal field theroy with central charge c — 1 (Ref. 3] 
and references therein). 

In the case of XXZ model a transverse field breaks the 



U(l) symmetry of the Hamiltonian to a lower, Ising- like, 
which develops a gap. The ground state then has long 
range anti- ferromagnetic order (0 < A < 1). However 
due to non-zero projection of order parameter on field 
axis it is a spin-flop Nee l state. In fact at a special field 
(hff = 2y / 2J x (J x + A)) the ground state is known ex- 
actly to be of classical Neel type Phase diagram, 
scaling of gap and some of the low excited states at /ijv 
has been studied in Ref. @ . The gap vanishes at a critical 
field h c , where a transition to paramagnetic phase occurs. 
Classical approach to this model reveals the mean field 
results 7 which is exact as s — ► 00. However the study 
of critical region needs quantum fluctuations to be taken 
into account. Exact diagonalization £ and Density Ma- 
trix Renormalization Group (DMRG) — gives us the prop- 
erties of stable phases. A bosonization approch to this 
model in certain limits leads to a nontrivial fixed point 
and a gapless line which separates two gapped phases^, 
moreover the connection to the axial next-nearest neigh- 
bor Ising model (ANNNI) has been addressed. The ap- 
plicability of mean-field approximation has been studied 
by comparing with the DMRG results of magnetization 
and structure factor 11 . Recently the effect of longitudinal 
magnetic field on both Ising model in Transverse Field 
(TF)ii and XXZ model in TF has been discussed^. Here 
we are going to present the phase diagram of XYZ model, 
Eq. JIJ, by means of RG flow of coupling constants to 
show explicitly its universality class. 

Apart from theoretical point of view, recent experi- 
ments on CS2C0CI4 in the presence of transverse mag- 
netic field can be explained by XYZ model with A = 0.25 
14 . Using Quantum Renormalization Group (QRG) we 
will show explicitly that the anisotropy is not relevant 
and the universality class is governed by Ising model in 
Transverse Field (ITF). In addition QRG results rule out 
the existence of spin liquid phase between spin flop and 
paramagnetic phases which are separated at the critical 
field (h c ). Exact diagonalization data supports our QRG 
results by calculating the structure factor and magneti- 
zation of finite chain sizes. Our results are in good agree- 
ment with the experimental data. We will also discuss 
on the reasons why magnetization does not saturate just 
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above critical point. 

Quantum RG scheme in real space is started by de- 
composing lattice into isolated blocks. The Hamiltonian 
of each block is diagonalized exactly and some of the low- 
lying states is kept to construct the basis for renormal- 
ized Hubert space. Finally the Hamiltonian is projected 
to the renormalized spaced. We have considered a two 
sites block and kept the ground ( | eO) ) and first (|el)) 
excited states of each block to construct the embedding 
operator (T = |el)(f | + |eO)(j. |) — . Energy eigenvalues 
are eO = - J x - J y - A and el = J x - ^4h 2 + (J y - A) 2 . 
The | t) and | J.) are renamed basis in the renormalized 
Hilbert space. The interaction between blocks define the 
effective interaction of renormalized chain where each 
block is considered as a single site. A remark is in order 
when projecting the Hamiltonian to the effective (renor- 
malized) Hilbert space. The effective Hamiltonian is not 
exactly similar to the initial one, i.e. the sign of cr^crV, 1 
and of a\ \ x terms is changed. To avoid this and pro- 
ducing a self similar Hamiltonian we first implement a 7r 
rotation around x-axis for even sites and leaves odd sites 
unchanged. Therefore the Hamiltonian is transformed to 
the following form, 

N/2 

H = . - Jytrlcrl+x ~ A« +1 - haf]. (2) 

i=l 

We note to interpret our final results in terms of this 
transformation. The renormalized Hamiltonian (H ren = 
T T H {trans f ormed)T) is similar to Eq.(0) with renor- 
malized coupling defined below. 

J x( {J y -A) 2 -d 2 ^ 2 

4 \{jy- A) 2 +l? 2 / 

= Jy (Jy-A + Of 

V 2 (jy-A) 2 +d 2 

eO-el J x s {J y - A) 2 -$ 2 \ 2 
2 2 \{J y - A) 2 +$ 2 ) 

A' A ^- A -^ 2 (3) 
" 2 {.J y -A) 2 +d 2 [6) 

where = ^4/j 2 + (J y - A) 2 - 2h. This RG-flow is not 
valid when h —> where the U(l) symmetry at J x = J y 
can not be recovered by Eq. J3J). It will be discussed 
later. However due to level crossing which happens for 
the eigenstates of block Hamiltonian, Eq. © is valid 
when g x < (1 + + 4gl)/2 and g A < 9x < !• This 
covers XYZ model (J x < J y ) in transverse field when 
< A < 1. The new parameters g x — 4^, gA — -r" an d 

Jy Jy 

gh = -T- are defined because these ratios actually define 

Jy 

competing phases. 

We have plotted the RG-flow (arrows) and different 
phases in Fig. ^ The RG equations (Eq. ((3J)) show 
running of A to zero. In other words the anisotropy 
term is irrelevant (0 < A < 1). So we have only plot- 
ted the A = plane, ft means that the universal- 
ity class of XYZ model in transverse field (TF) is the 




FIG. 1: Phase diagram of XY model in transverse field. Ar- 
rows show running of couplings under RG. Filled circles show 
fixed points and the open circle is the tricritical point where 
lines separating different phases merge. The tick line at the 
top of phase diagram (g x = +oo) is a line of fixed pointa—. 
Phase (I) is antiferromagnetic Ising in y-direction (spin-flop), 
(II) paramagnetic in x-direction and (III) is antiferromagnet 
in x-direction. is the ITF fixed point. 



same as XY model in TF. Moreover the exchange inter- 
action in the x-direction is also irrelevant while J x < J y . 
As J x vanishes under RG, there are only two effective 
terms in the Hamiltonian. This is exactly the case of 
Ising model in TF (ITF model). So the interplay of 
J v <jfaf +1 and h(af + of +1 ) defines either ordering in y 
or paramagnetic in x direction. Solving the RG equa- 
tion for fixed points, we found the non-trivial fixed point 
g* t = (g x = 0,<7/j 1.26,5a = 0) apart from the other 
which is &t(g x = 0,<?/j = oo,<?a = 0) and represents sat- 
urated ferromagnet. We have linearised the RG-flow at 
g^ and found one relevant direction (whose eigenvalue 
is larger than one). The eigenvalues and correspond- 
ing eigenvectors of linearised RG at <? ; * in [g x ,gh,gi\) 
space are: |A X = 1.59) = (0, 1, 0); |A 2 = 0.31) = 
(1, 1.64, 0); |A 3 = 0.46) = (0, 0.62, 1). The relevant di- 
rection ( | Ai ) ) is the horizontal line passing through g* h 
and | A2) is the tick line ending at g^. The critical ex- 
ponents at this fixed point are (3 = 0.41,^ = 1.48 and 
z = 0.55. The discrepancies of exponents from exact val- 
ues (/? = 0.125, = 1 and z = l)i£ are the result of 
2-sites blocking, however these are exactly equal to the 
exponents of ITF chain which is calculated by QRG 16 . 
As far as g x < 1, the control parameter is gh- When 
9h < 9h (phase (I)), the staggered magnetization in y- 
direction (SM y ) is nonzero which is the order parameter 
to represents the phase transition at g^ (the line which 
ends at g^). However magnetization in x-direction (M x ) 
is also nonzero which causes to consider this phase as a 
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FIG. 2: a The order parameter (SM y ) and magnetization in 
x-direction (M x ) versus transverse field. QRG and extrap- 
olated Lanczos results are compared for (M x ). b Lanczos 
results of M x vs. transverse field for N = 12, 16, 20, 24 and 
extrapolation to N — » oo, at A = 0.25 and J x = J y = 1. 
The inset shows that M x behaves qualitatively the same for 
A = and 0.25. 



spin flop phase. This is an Ising like phase which has a 
nonzero gap. This gap is going to be closed at g^ where 
the transition to paramagnetic phase takes place. At this 
point the quantum fluctuation of TF destroys the antifer- 
romagnetic (AF) ordering completely. The paramagnetic 
phase (II) appears at gh > g^ where spins are aligned in 
the field direction and will be saturated in high TF. Note 
that the proper order parameter for this phase transition 
is staggered magnetization in y-direction. So it is not nec- 
essary to gain the saturation value for M x just after g^. 
This also happens in ITF model. We have plotted both 
SM y and M x in Fig.^a). The comparison with Lanczos 
results show very good qualitative agreement. Although 
it is not expected that QRG gives good quantitative re- 
sults we got fairly well agreement with Lanczos results. 

To discuss the behavior close to h = 0, we need to take 
into account the U(l) symmetry in the QRG scheme. So 
we will consider the XY model at h = and the effect 
of TF is taken into account by perturbation. In this case 
the only relevant parameter is g x . Implementing a 3 sites 
blocking, the RG flow is: g' x = g x , which has 2 stable 
g* = 0,oo and an unstable fixed point g* = 1. The 
stable fixed points define two AF Ising phases ordered 
in y-direction (g* = 0) and x-direction (g* — oo). The 



g* = 1 is the critical point where a transition occurs 
between two stable phases. Now the transverse field is 
considered pcrturbatively which gives the following RG 
flow for gh- 

, / 2gsyT+ gf- g| \ n ,.s 

9h = { i + tf ) 9h ; 9h ^° (4) 

The perturbation approach is justified since gh — > 0. For 
any value of g x , Eq. leads to g' h < gh, which means 
the direction of flow is toward the g x axis. As a result of 
QRG at gh — we expect to have a phase transition at 
small gh by changing g x close to g x ~ 1. The boundary 
of this phase transition is shown by dashed line in Fig.^ 
which is the gapless line reported in Ref . [Iol| . This line 
represents the phase transition between phases (I) and 
(III), AF Ising in y- and x-direction, respectively. As 
g x — > oo (J y — > 0) the model behaves as an AF Ising in 
a longitudinal magnetic field. In this limit a first order 
phase transition at -j- = 1 divides the AF -j- < 1 from 

paramagnetic -j- > 1 phases. A line of fixed points comes 
out of a 3-sites block QRG 18 for j- < 1 which has been 
shown as a tick line at the top of phase diagram (Fig.Q. 
Thus a line with slope |j = 1 (as J y — > 0) constructs 
the boundary of phase transition between (II) and (III). 
This phase transition is in the universality class of AF 
Ising in a magnetic field. To complete the structure of 
phase diagram we propose a tri-critical point (open circle 
in Fig. which is the coexistence point of three phases. 
Still we do not have an RG equation at this point. 

We have implemented the Lanczos algorithm on finite 
sizes (N — 12, 16, 20, 24) using periodic boundary con- 
ditions to calculate M x and structure factors both in x 
and y directions. In Fig. \^b) we have plotted M x for 
different chain sizes at A = 0.25 and an extrapolation to 
N — * oo. The value of A = 0.25 is chosen to fit the case 
of CS2C0CI4. The general behavior is similar to what we 
have obtained from QRG (Fig.|^a)). There is no sharp 
transition to the saturation value at a given h because 
M x is not the proper order parameter to this phase tran- 
sition. Oscillations of M x at finite N for h < h c is the 
result of level crossing between ground and first excited 
states of this model. The last level crossing happens at 
h . We have also plotted the case of A = to show 
the same qualitative behavior as A = 0.25 in the inset 
of Fig. [2Kb) . Lanczos results leads to SM y — for any 
value of h, since in a finite system no symmetry break- 
ing happens. However the structure factor (S yv (q = ir)) 
diverges in the ordered phase as N — ► 00. The structure 
factor at momentum q is defined as 

S aa (q) = Y,<°S< >e iqr ; oc = x, y (5) 

r 

In Fig. Eta), S vv (q = tt) is plotted versus for different 
transverse field. As far as h > 3.1, S vv (q = tt) grows 
slowly and shows saturation at a finite value when N — > 
00. In the other hand a super linear behavior versus N 
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FIG. 3: Structure factor (a) S yy (q = n), (b) S xx (q = 0) 
versus N for different transverse field. S vv (q — 7r) shows di- 
vergence as N — » oo while h < h c — 3.1 (in the ordered 
phase). All plots for S xx (q = 0) show divergence in thermo- 
dynamic limit (N — > oo). However super linear behavior for 
h < h c — 3.1 and almost linear behavior for h > h a is the sign 
of two different phases. 



shows a divergence of structure factor for h < 3.1. It 
corresponds to ordered phase which is AF in y-direction. 
Thus the critical field at A = 0.25 is h c = 3.1 ± 0.05. A 
similar computation results to h c = 2.9 ± 0.05 for A = 
0. To get an impression that the QRG results are very 
surprising we just mention the value of critical field for 
comparison with Lanczos ones, h c (A = 0.25) = 3.32 and 
h c (A = 0) = 3.12. 

We have also plotted the structure factor S xx {q = 0) 
versus N in Fig. Efb). This shows divergence for any 



value of h as N — > oo which verifies ordering in x- 
direction. The spin flop phase (I) has nonzero M x which 
increases by h to the saturation value in paramagnetic 
phase (II). However we observe different qualitative be- 
havior for h < h c = 3.1 and h > h c . The former is super 
linear and the latter is almost linear. As mentioned be- 
fore, M x is not the proper order parameter and is not 
expected to be saturated at a specific h. The saturation 
happens for enough large value of TF. 

Summing up the QRG and numerical results, we claim 
that the universality class of XYZ model in TF (0 < A < 
1) is the ITF model. Thus there exists only two stable 
phases, namely (I) and (II), which are distinguished by a 
critical field at h c . In this respect there is no spin liquid 
phase just after transition point. We found very good 
agreement in the sense of universal behavior with the ex- 
perimental results — on Cs-iCoGl^. We have obtained 
the corresponding critical magnetic field H c = 1.3 T com- 
paring with the reported H c = 2.1 T . The difference 
should come from two doublets nature (s = 3/2) of ac- 
tual material and the effective Hamiltonian of s = 1/2 in 
our calculation which is responsible for low fields. The 
other mismatching is the observed crossover behavior in 
M x . As proposed in Ref— the crossover behavior is re- 
lated to the saturation of the lower doublet of Co 2+ and 
the inset of higher doublet effects. However for the XYZ 
chain as a spin 1/2 model this does not happen. At 
J x = J y , applying small noncommuting fields break the 
17(1) rotational symmetry and develops a gap which has 
the consequence of promoting long-range order in a spin- 
flop phase (I). Increasing field stabilizes the perpendicu- 
lar AF order which can be observed by the maximum in 
SM y . Higher TF reduces ordering up to a critical field 
(h c ) where gap vanishes. Just after this transition point 
a gapped paramagnetic phase appears (II). 
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